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The derivation of the Feynman path integral based on the Trotter product formula 
is extended to the case where the system is in a magnetic field. 



^ ■ I. INTRODUCTION 

qh| There are many ways to derive the path integral. In his original paper, Feynmani gave 

a physical motivation for a particular form and verified its correctness by propagating the 
wave function for an infinitesimal time step. In Ref. 0the opposite perspective is used: start 
with a Hamiltonian and use the Trotter product formula^ to arrive at the sum-over-paths. 
\ Some of Feynman's poetry is lost, but one connects more easily to familiar formulations. 
The Trotter-product-formula-based derivation in Ref. is limited to Hamiltonians of the 
form H = /2m + V{x), and skirts the more delicate case of a particle in a magnetic field. 
As is well known, the latter situation requires that when writing the classical action for the 
broken line path in the path integral, the vector potential, A{x), be approximated by its 
average at the ends of each straight segment of the path. One can also evaluate A at the 
midpoint of that segment, giving rise to the "midpoint rule." This feature is one of the 
fundamental properties of the path integral, reflecting the close relation between the path 
integral and the Wiener integral, and the "infinite velocities" or non-differentiability of the 
paths in both cases. Yet more delicate evaluations are necessary for the path integral on 
curved spaces, so it would be useful to understand how to deal with these rough paths from 
the Trotter derivation perspective. 

In this article we use a slight variation of the derivation in Ref. to show how we can 
arrive at the correct path integral, including a magnetic field. Some years ago two of us 
published such a derivation,— but it was phrased in the language of stochastic processes. 
The present work uses a framework more familiar to physicists and is more streamlined. 
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The earlier paper includes path integrals on curved spaces, which are not treated here. 
Because the purpose of this article is to extend familiar path integral techniques, we do not 
elaborate on how the formulae developed here make their way into the usual path integral 
representations. 

II. BREAKING UP THE PROPAGATOR 

The propagator is the operator exp{—iHt/h), with H the Hamiltonian. H is taken to be 

H=^{p--^A{x)r + V{x), (1) 

in the usual notation. For a, 6 G M'^, the propagator is exactly given by 

G{b,t;a) = {b\exp{-iHt/h)\a) (2) 
= {b\[exp{-tHt/Nh)]^\a) (3) 

= {b\exp{-iHt/Nh) j d'^Xi\xi){xi\exp{-iHt/Nh) 
X j d^X2\x2){x2\exp{—iHt/Nh) X . . . 

X J d^XN-i\x]\f-i){xN-i\exp{-iHt/Nh)\a). (4) 
Equation ^ can be written more concisely as 

Yl d^Xk Yl Gixe+i, e; Xe), (5) 

k=l 1=0 

with e = t/N, X]\f = b, and Xq = a. 

The basic fact of life about the short time propagators, G(xe+i, e; xi), is that for N oo, 
we need only maintain 0(e) accuracy. That is, they can be replaced by other objects that 
differ only by terms that go to zero faster than e. For ordinary numbers this can be seen 
by recalling that e^' = limiv-^oo(l + f/)'^ y by observing that adding to x any quantity 
that goes to zero with does not change the result. That is, (1 + ^'"^^ ) ^ provided 
aj^ —>■ 0. To show that this limit works for unbounded operators is of course a much bigger 
project and the reader is referred to Refs. iandi and to references therein. 

Our goal therefore is to approximate G{x, e; y) to first order in e. We return to operator 
language. The starting point is to observe that for two operators 

\2 

exp[X{A + B)] =exp{\A)exp{XB)exp{ — [B,A]+0{\^)). (6) 

The propagator is of this form, with X = e, A = —iK/h, B = ~iV/h, and K the kinetic 
energy. We can also write Eq. (0) as^ 

2^ 



exp[A(A + B)] = exp(AA) exp(A5) + O(A'). (7) 



This formula lets us replace exp[A(y4 + B)] by exp(Ay4) exp(Ai?) when A is e. 
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In the absence of a vector potential one term is a function of momentum only and the 
other a function of position only, allowing the breakup that leads to the usual path integral. 
For later reference we review this procedure. The potential, V, is diagonal in position space, 
so that we immediately obtain 

G{x,e;y) = {x\exp{~iKe/h)\y) exp{-iV{y)e/h) + 0{e^). (8) 

If K is simply p^/2m, the next step is straightforward. By inserting J d^p \p){p\ to the left 
of \y), we find that exp{—iKe/h) gives- 

|d'p{x\exp{-^p\/h2m)\p){p\y) = {^f'exp (^y^^^)- (9) 

(See Eq. (|T7|) for the form of the momentum eigenstates.) If we combine Eqs. (jS)) and Q, 
we obtain the expression, which when iterated, gives the path integral in three dimensions: 

^2 



"^^--■'y^-^^E^r^^piiC^^-'Viy^ H-o(.=). (10) 

Thus, if we insert Eq. fllO|) in Eq. (jSJ, we obtain the classical action in the exponent. Note 
that by interchanging K and V in Eq. ((Tj), the argument of V in Eq. pO|) becomes x rather 
than y. As indicated, this changes the short-time propagator by less than 0(e) and therefore 
does not change the final result. Indeed, any appropriately weighted sum of V{x) and V{y) 
can be used, or as we shall see below, V can be taken at any point on the line between x 
and y. 



I ,'m[x — y 



III. UNCOMPLETING THE SQUARE 

When a magnetic field is present, the expression exp{—iKe/h) is more complicated, be- 
cause K = [p — eA{x) / /2m. Under these circumstances, inserting a momentum- state 
resolution of the identity is inadequate, because A is a function of x. The idea of the next 
step is to look at the square root of , which because it is a sum of p and A, can be resolved 
by separately expressing them in momentum and position space bases. Taking the square 
root in the exponent is accomplished at the expense of introducing an additional integral, a 
process variously known as "uncompleting the square'^ or the "Gaussian trick"—. The idea 
is based on the identity 

exp (^) = (-1=) Y d'uex^i-'^ + h-u). (11) 

We use the following variant of Eq. : 

exp ( — ie-^—) = ( }. Y [ d^u exp (i— i\ b ■ u). (12) 

For convenience, we let a{x) = eA{x)/c. We want to evaluate Eq. (jH)) and concentrate on 
the kinetic energy part, which, with the aid of Eq. (fT^ is rewritten as 

GK{x,e;y) = {x\exp{-iKe/h)\y) 

Xp ,y 

\mh 

a)-u)\y). (13) 
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Equation ()13p shows where the trouble enters. The p has been separated from the a with 
which it does not commute. That's the good part. Nevertheless, Eq. ((Tj) would not provide 
the needed accuracy, because it is the square root of e that appears in Eq. (fT^ (A of Eq. ((Tj) 
is now ^/e). The error is the square of this quantity, namely e itself, which cannot be 
neglected. This problem will occur whenever an expression mixes non-commuting variables, 
for example in a curved-space metric (or kinetic energy). Generally speaking, this is why 
the operator ordering problem does not get solved in the path integral. The way to deal 
with this problem is to improve on Eq. ((Tj). To this we now turn. 



IV. A MORE ACCURATE PRODUCT FORMULA 



A slight variation of Eq. ((Tj) provides 0(A'^) accuracy: 

exp[X{A + B)] = exp{XB/2) exp(AA) exp(A5/2) + 0{X^). 



(14) 



Equation (fT^ can be checked by direct expansion of the exponents-^ A different symmetriza- 
tion was used in Ref. 0- Equation (fT^ has also been used for increased accuracy in the 
numerical evaluation of path integrals.— 

We apply Eq. (fT^ to (a:;| exp ( — i\fe{p — a) ■ u/\fmh)j \y) to obtain 



{x\exp[-iJ — {p-a) ■u)\y) 



{x\ exp (zW — -a ■ u/i) 

^ mn ' 



X exp ( — i\l — -p ■ u) exp (i^l — -a ■ u/2) \y) 
^ mn v mn 



exp (i\ — ^o-ix) ■ u/2) 
^ mn ' 



x(a;|exp ( - i^l-^p-u)\y) exp {iJ—a{y) ■ u/2). 



(15) 



The error in the first "equality" is 0{e^^'^), although it is not explicitly indicated. Once again, 
a momentum resolution of the identity is inserted in the p-dependent portion of Eq. ()15j): 

i.l.M-^^r ") I.) ^/'^V (.1 exp i-^^P-u) (16) 
If we recall that 

{x\p) = {2nh)~^^^ exp{ip ■ x/h), (17) 

then Eq. (fT^ becomes 

{x\ exp(— ii/ep ■ u / V mh)\y) = (--^^ [ d^p exp{—i\^p ■ u/y/mh + ip ■ {x — y)/h) 

VzTrn,/ J 



1 



6^{-y/eu/Vmh + {x - y)/h) 



mh 



6^{—u + {x — y)\/m/he) 



(18) 



5 



— u 

m 



The integral over u in Eq. ()13|) is now trivial, yielding for Gk-, 

with a = eA/c. For the classical Lagrangian, a magnetic field contributes v ■ eA{x) / c. If we 
multiply and divide by e in the last expression in Eq. (|T9|l . we obtain exactly the appropriate 
term. For completeness, we restore A and V , yielding, 

m ii fz ^m{x-yY [x - y) e,A{x) + A{y) 



+ 0(e3/2) ^ (^20) 

The argument of the exponent is seen to be ie/h times the classical Lagrangian. 

Note that what naturally arises is the average of the vector potential. A, at the endpoints 
of the broken line path. The "midpoint" way of attaining the same level of precision uses 
A [(a; + y)/2]. The difference between these forms is essentially a second derivative of A 
times {x — y)"^. The latter is of order e and in turn multiplies an additional power of (x — y), 
so that the difference, of order e^/^, can be neglected. (This observation also lies behind our 
remark in Sec. Othat V can be evaluated anywhere along the line between x and y. In 
Eq. (|2Up a similar freedom exists for V^s argument.) 



V. DISCUSSION 

The beauty and efficacy of the path integral come at a price: ordinary calculus can no 
longer be taken for granted. As for the Wiener integral used for Brownian motion, the 
"paths" one sums over are rough. They are nowhere differentiable and if one allows time 
differences on a trajectory to go to zero, the corresponding spatial differences go to zero far 
more slowly, to wit. Ax ~ \/At. For this reason, when evaluating the action on broken 
line paths when a magnetic field is present, the term Ax ■ A must be handled carefully. 
Feynmani introduced the midpoint rule, namely the use of {x — y) ■ A[[x + y)/2]. 

In this article we have shown how this prescription arises from a derivation of the path 
integral that begins from the Hamiltonian and uses a slightly more accurate form of the 
Trotter product formula. Similar arguments can be made for path integrals on curved 
spaces,— where the need for accuracy is even greater.— 

We emphasize that this feature of the path integral is not some minor annoyance, an 
inessential complication. It is a central feature of quantum mechanics, evident for example 
in the fact that velocity cannot be defined for quantum systems (at least not as the limit 
of Ax /At). Moreover, it plays a central role in applications of the Wiener integral, for 
example, in the derivation of the Black-Scholes formula.— 
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